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\^ ' Abstract 

o '. 

tH- , We analyze the extraction of weak phases from CP violation in B ^> 

O I TT~^TT~ decays. We propose to determine the unitarity triangle {p, ff) by 

combining the information on mixing induced CP violation in B -^ vr+vr^, 
Mh, S, with the precision observable sin 2/3 obtained from the CP asymmetry 

Q^l in i? ^ ipKs- It is then possible to write down exact analytical expressions 

^ ' for p and f) as simple functions of the observables S and sin 2(3, and of the 

penguin parameters r and (j). As an application clean lower bounds on f] 
and 1 — /? can be derived as functions of S and sin 2/3, essentially without 
^ ' hadronic uncertainty. Computing r and (j) within QCD factorization yields 

j^ ■ precise determinations of p and f/ since the dependence on r and (f) is rather 

weak. It is emphasized that the sensitivity to the phase (p enters only at 
second order and is extremely small for moderate values of this phase, 
predicted in the heavy-quark limit. Transparent analytical formulas are 
further given and discussed for the parameter C of direct CP violation 
in i? — > ■n~^'K~. We also discuss alternative ways to analyze S and C 
that can be useful if new physics affects B^-B^i mixing. Predictions and 
uncertainties for r and (p in QCD factorization are examined in detail. It 
is pointed out that a simultaneous expansion in l/rrif, and 1/N leads to 
interesting simplifications. At first order infrared divergences are absent, 
while the most important effects are retained. Independent experimental 
tests of the factorization framework are briefly discussed. 
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1 Introduction 

The main goal of the current experimental program at the SLAC and KEK B- 
meson factories is a stringent test of the standard model description of CP vio- 
lation. In the future this aim will be pursued with measurements of still higher 
precision from hadron machines at Fermilab and CERN. A crucial benchmark 
is the time- dependent CP violation in i? — » ip^s decays, which allows us to in- 
fer the CKM phase f3 with negligible hadronic uncertainties. Likewise of central 
importance for obtaining additional information on CKM parameters is the time- 
dependent CP violation, both mixing-induced (S) and direct (C), in B ^ tt+tt". 
However, in this case the extraction of weak phases is complicated by a penguin 
component in the decay amplitude, which carries a weak phase different from the 
leading tree-level contribution. This leads to a dependence of the CP asymmetries 
in B ^ 7T^TT~ on hadronic physics and to corresponding theoretical uncertainties. 
A possible strategy to circumvent this problem is the isospin analysis by Gronau 
and London JJ, where also the branching ratios of B~^ -^ tt+tt", B -^ 7r°7r° and 
their charge conjugates have to be measured. This method is theoretically very 
clean, but the difficulty to measure B -^ ti^tt^ decays with sufficient accuracy and 
the existence of discrete ambiguities are likely to prevent a successful realization. 
It is the purpose of this paper to demonstrate how the information on weak 
phases contained in the CP asymmetries of i? ^ tt^tt^ itself can be extracted in 
an optimal way. To some extent theoretical input on the penguin-to-tree ratio 
will be needed and can be provided by the QCD factorization approach 



However, we will show that the impact of uncertainties in the calculation is in 
fact very mild. Moreover, even if the detailed predictions of QCD factorization 
are ignored, it is still possible to derive rigorous bounds on the CKM unitarity 
triangle, using only very conservative assumptions. 

In order to derive these results we propose the following strategy. First, the 
time-dependent CP asymmetries in B ^ vr+vr" are expressed in terms of Wolfen- 
stein parameters p and f]. At the same time purely hadronic quantities are sys- 
tematically isolated from CKM parameters, introducing magnitude r and phase 
of a suitably normalized penguin-to-tree-ratio |3]. We then combine the observ- 
able S{p, fj, r, (f)) with the accurately known value of sin 2/3(p, fj) from B -^ ipKs- 
This allows us to obtain the exact unitarity triangle, p and f/, in a simple ana- 
lytical form, depending only on sin 2/5, S and the hadronic quantities r and (j). 
The dependence on the latter turns out to be particularly transparent, which 
greatly facilitates any further analysis. We are then able to derive bounds on the 
unitarity triangle practically free of hadronic uncertainties, or to fix p and f/ with 
theoretical input for r and 0. 

There is already an extensive literature on the subject of extracting infor- 
mation on weak mixing angles from CP violation in 5 — > 7r"'"7r~ |HMl6j. In these 
papers important aspects of the problem have been discussed and suggestion were 
made to constrain theoretical uncertainties. Here we present a new way of ex- 
ploiting the information contained in the CP violation observables S and sin 2(3. 



The crucial elements are a definition of hadronic quantities r and independent 
of CKM parameters, the direct formulation of weak phases in terms of the basic 
Wolfenstein parameters p, f/, the resulting analytical determination of the unitar- 
ity triangle and the exact, explicit and very simple dependence on r and 0. This 
in turn greatly facilitates the analysis of theoretical uncertainties and gives, in 
combination with results based on the heavy-quark limit, robust determinations 
of the unitarity triangle, or rigorous CKM bounds with minimal hadronic input. 
These ideas were first presented in [T7]. Subsequently, this analysis has been 
further discussed by Botella and Silva [iHj and Lavoura [T9] . 

The paper is organized as follows. In sec. 2 we collect important basic formu- 
las describing CP violation in i? — > vr+vr". In sec. 3 we discuss the theory of the 
penguin parameters r and in the framework of QCD factorization. Based on 
previous work we address in particular the issue of theoretical uncertainties. In 
addition we investigate the analysis oi B ^ tt+tt^ decay amplitudes in a simul- 
taneous expansion in both l/rrih and 1/A^, where A^ is the number of colours. An 
interesting pattern of systematic simplifications resulting from the double expan- 
sion is pointed out. After this discussion of the hadronic input, we turn to our 
phenomenological analysis. Sec. 4 explores the determination of the unitarity 
triangle from S and sin 2/5 within the standard model. Simple analytical expres- 
sions are presented and theoretically clean bounds are derived. We also evaluate 
the standard model expectation for S using results from QCD factorization. Sect. 
5 considers the determination of CKM angle 7 from S, sin 2/5 and iKib/Kbl. This 
alternative possibility is useful if new physics affects the phase of B^-Bd mixing. 
Sec. 6 examines what can be learned from C, the parameter of direct CP viola- 
tion in S — > vr+TT". Methods to validate the predictions of QCD factorization for 
B -^ TT+TT^ using additional observables are described in sec. 7. We summarize 
our main results in sec. 8. 



2 Basic Formulas 

The time-dependent CP asymmetry in B -^ tt+tt^ decays is defined by 

B{B{t) -^ TT + TT-) - B{B{t) -^ TT + TT-) 



^cpW 



B{B{t) -^ TT+TT-) + B{B{t) -^ TT + TT-) 

-S sin(Amijt) + C cos(AmBt) (1) 



where 



1 + leP 1 + leP e+'^ + P/T ^ ' 

In terms of the Wolfenstein parameters p and f/ |2UI21j the CKM phase factors 

read 

±i^ ^ P^iy _2i^ ^ (1 - p)^ -ff - 2if]{l - p) 



The penguin-to-tree ratio P/T can be written as 



T Vp^ + V 



7<1 ^^2 



(4) 



The real parameters r and (f) defined in this way are pure strong interaction 
quantities without further dependence on CKM variables. 

For any given values of r and a measurement of S defines a curve in the (p, 
f/)-plane. Using the relations above this constraint is given by the equation 

_ 2f][p^ + ff - P - p(l -P) + (p^ + ff - l)rcos0] 
{{1 — pY + ff){p^ + f]'^ + r'^ + 2rpcos(j)) 

Similarly the relation between C and p, fj reads 

C = '^^^ ^^^ ^ (6) 

p2 _|_ ^2 _|_ j,2 _|_ 2rpcos(j) 

This is equivalent to 

r sin (j)\ /I 



(p + rcos0)^+ U-^^^ 1 = (^ -Ij (rsin0)^ (7) 

describing a circle in the (p, f/)-plane with centre at (— rcos0, (rsin0)/C) and 
radius r sin (pyjl — C^ jC . 

The current experimental results for S and C are 

^ = +0.02 ±0.34 ±0.05 (BaBar|22I) -1.00 ± 0.21 ± 0.07 (Belle[2SI) 

(8) 
C = -0.30 ±0.25 ±0.04 (BaBar[22]) -0.58 ± 0.15 ± 0.07 (Belle[2ni) 

A recent preliminary update from BaBar gives |21] 

5 = -0.40 ±0.22 ±0.03 C = -0.19 ±0.19 ±0.05 (9) 

3 Penguin Contribution 

In this section we discuss theoretical calculations of the penguin contribution in 
B -^ TT+TT^. The analysis is based on the effective weak hamiltonian 

^e// = %Y.\ ( Ci g? ± C72 gf ± E a g. I ± fix. (lo) 

V2p=u,c \ i=3,...,6,8g J 

where Ci are Wilson coefficients, known at next-to-leading order |25j, and Qi 
{i = 1, . . . , 6) are local four-quark operators with fiavour structure {db){qq), q = 
u, d, s, c, b. Qsg is the chromomagnetic operator ~ nibda ■ G(l ± 75)6. The CKM 
factors are here denoted by Xp = V*^Vpb. 



3.1 QCD Factorization 

The penguin parameter r e**^ has been computed in Pj in the framework of QCD 
factorization. The result can be expressed in the form 

ai + al + rlal + rA[hi + h + 2hi] 

where we neglected the very small effects from electroweak penguin operators. 
The factorization coefficients a, are linear combinations of the Wilson coefficients 
Ci in the effective weak Hamiltonian and include the 0{as) corrections from hard 
gluon interactions in the weak matrix elements. Their expressions can be found 
in ^. The quantities r^ and ta are defined by 



2m^ fsU 

ty A 



rW = ,^....,^ ...:, ,^...^^ rA = , . Z'Z.... (12) 



r^ is defined in terms of the MS quark masses mq{fi) and depends on the renor- 
malization scale /x. -Fq^^'^(O) is a -B ^ tt transition form factor, evaluated at 
momentum transfer g^ = m^ ~ 0. 

Both quantities in (fT^ are of subleading power r'!^ ^ va ^ KqcDl'fnb- ^a ~ 
0.003 is numerically very small. It sets the scale for the weak annihilation effects 
in the amplitude, which are parametrized by the hi in (jlip [4 . They represent 
power corrections that are not calculable in QCD factorization. Model-dependent 
estimates for these subleading effects have been given in [^ in order to assess the 
corresponding uncertainties. On the other hand, r^(1.5 GeV) ~ 0.7 is numerically 
sizable. Still the important penguin contributions ag, p = u, c, are calculable and 
can be included in the analysis. A third class of power corrections that need to 
be considered are uncalculable spectator interactions, some of which also come 
with the parameter r^. These enter the a^ in (fTT|) and were also estimated in [3j. 

In Tabled we show the values for r and from a calculation within the QCD 
factorization framework as described in |1]. We also display the uncertainties 
from various sources, distinguishing two classes. In the upper part we give the 
uncertainties from input into the factorization formulas at next-to-leading order, 
as well as the sensitivity to the renormalization scale /i. This input is defined 
in Table |21 The second class of uncertainty is due to the model estimates 
employed for power corrections. As in |1] these effects are parameterized by 
phenomenological quantities 

XH,A=(l+PH,Ae"^"'')i^^ (13) 

^ ^ Aft 

that enter power corrections to hard spectator scattering (if) and weak annihila- 
tion effects (A). The default values have ph,a = 0. They depend on an infra-red 
cut-off parameter Ah, which we take as A^ = 0.5 GeV. An error of 100% is then 
assigned to this estimate by allowing for arbitrary phases (pn, (pA and taking pn, 
Pa between and 1. The impact on r and (p of this second class of uncertainties 





/i 


rriu + md 


rric 


fB 


^0 


«2 


Ab 


r = 0.107 


±0.005 


±0.019 


±0.002 


±0.003 


±0.002 


±0.003 


±0.002 


= 0.150 


±0.023 


±0.001 


±0.057 


±0.002 


±0.002 


±0.003 


±0.002 





(PhAh) 


(pa, 0a) 


r = 0.107 


±0.001 


±0.024 


= 0.150 


±0.010 


±0.24 



Table 1: Theoretical values for r and and their uncertainties from various 
sources within QCD factorization. The upper part displays uncertainties from 
input into the factorization formulas. The lower part gives the uncertainties 
from a model estimate of power corrections (see text for details). 



rriu + md 


m^imb) 


fB 


Fo^--(O) 


«2 


Ab 


0.0091 ±0.0021 


1.3 ±0.2 


0.18 ±0.04 


0.28 ±0.05 


0.1 ±0.3 


0.35 ±0.15 



Table 2: Input used for Table U\ We take /i G [mb/2, 2mh] . The values for 
^u ± ^TT'd = iji^u ± "^d)(2GeV), rricimi,), fs and Xb are in GeV. 

is shown in the lower part of Table C] and is seen to be completely dominated by 
the annihilation contributions. 

Adding the errors in quadrature we find 



0.107 ±0.020 ±0.024 
0.15 ±0.06 ±0.24 



(14) 
(15) 



where the first (second) errors are from the first (second) class of uncertainties. 
Combining both in quadrature we finally arrive at 

r = 0.107 ±0.031 = 0.15 ±0.25 (16) 

which we take as our reference predictions for r and in QCD factorization. 



3.2 Expansion in l/rrib and 1/A^ 

In order to obtain additional insight into the structure of hadronic i?-decay ampli- 
tudes, it will be interesting to consider these quantities in a simultaneous expan- 
sion in powers of 1/mi, and 1/A^, where A^ is the number of colours. Expanding in 
l/rrib alone corresponds to the framework of QCD factorization, implying naive 
factorization at leading order, which receives calculable corrections. Large- A^ ex- 
pansion of weak decay amplitudes gives an entirely different justification for naive 



factorization, which holds at leading order in 1/A^. The large-A^ limit has previ- 
ously been applied to weak decays of kaons and D mesons [26 27 28, . Here we 
would like to explore the consequences of combining the heavy-quark and large- iV 
limits in the analysis of important subleading corrections to naive factorization. 

For this purpose we treat the Wilson coefficients Ci, C2, C4, Cg and C^g as 
quantities of order one. In the limit N ^ 00 strictly speaking only Ci and Cgg 
are nonvanishing. However, C2, C4 and Cq vanish slower than 1/N, a.s N —* 00, 
if the large logarithm InMw/mfy is considered ~ N, in accordance with the usual 
renormalization group (RG) counting aglTiMw/iTib ~ 1 and with as ~ 1/A^. 
More specifically, C2, C^ ~ hiN/N and Cq ~ 1/N'^/^'^. The formal treatment of 
these coefficients as order unity is identical to the usual counting of the coefficients 
in RG improved perturbation theory. From this latter case it is also clear that 
the small numerical size of C2 and especially of the penguin coefficients C4, Cq 
is related to small anomalous dimensions, which are small accidentally, but not 
because of a particular parametric suppression. The coefficients C3, C5 on the 
other hand are suppressed relative to C4 and Cg by an explicit factor of N . 
We will thus take C3, C5 ~ 1/A^. Similar considerations can be found in J2Z1 
1^ . Concerning the heavy-quark limit uii, ^ ^qcd there is no difference to the 
conventional counting in RG improved perturbation theory, where Mw ^ fni, is 
assumed. 

After discussing the Wilson coefficients we next turn to the hadronic matrix 
elements in QCD factorization. We first need to determine how the various quan- 
tities entering these matrix elements scale for large rrih and A^. Based on these 
results we shall then expand the factorization coefficients Oj and 6j to next-to- 
leading order in a double expansion in l/nib and 1/N. That is, we keep terms of 
order one, as well as corrections suppressed by either a single power of l/rrif, or 
1/A^. We neglect terms that exhibit a suppression by two or more powers of the 
expansion parameters l/rrib or 1/A^, such as 1/ml, l/N"^ and l/rribN. It turns 
out that in this approximation all subleading contributions suffering from in- 
frared endpoint singularities in the QCD factorization approach are absent. This 
includes both spectator interactions of subleading twist and all weak annihilation 
amplitudes. On the other hand, important effects as hard QCD corrections or the 
penguin contribution from Qq, which is formally power suppressed, are retained. 

Let us postpone annihilation effects for the moment and examine first the 
coefficients a,. The most important ones for our purpose are ai, 04 and Og, which 
may be written as |lj 



ai 


= 


^ .C2 


al 


= 


'^-N 


< 


= 


'^-N 


K, PI2 


PI 


3 are calc 



1 + 



1 + 



Cpa, 

Air 
Cpoii 



K 



+ 



6 



47r 
Cpd 



V. 



C2 CpT^Ois 

iV N 

CpOis PiT,2 



Atx 



+ 



47r 
Cpoi 



H„ 



C3 CpTTOis jj 
+ -^^ TT SiiT 



N N N 



pp 

s -r 7r,3 



An N 



(17) 



(19) 



are calculable quantities of order one (in both l/rub and 1/N), 
containing convolution integrals over pion light-cone distribution amplitudes ^. 



The coefficient iif,r7r describes hard spectator scattering and reads 



Ht, 



fsflT 



di 

m|Fo^--(0)7o i 



MO 



f' 


dx , 

X 


.{x) 


r 


dy 


1 




L 


y 



X 



My) + rZ-My) 



X 



(20) 



Here (pB is the leading-twist hght-cone distribution amphtude of the B meson, 
07r the one of the pion. 4>p{y) = 1 is the two-particle, twist-3 component of the 
pion light-cone wave function. We recall that the correction ~ ^Z4>p (defined as 
/f^7r,3 in (EOl)) is uncalculable, as indicated by the end-point divergence in the 
y-integral, but it is power suppressed in l/rrif,. 

The large-rrib, large- A^ scaling of the various terms is as follows: 



^B, 



We then have 



Further we note Cp 
and 1/A^ we find 



^'^(O) 



~ 1 
1/m 



^~l/mfe x,y 



3/2 
b 



fB ~ N'/y 



m 



1/2 



/. - iVV^ 



-f^7r7r,2 ~ 



N 



H. 



TTTT.S 



N/rUb 



(21) 
(22) 

(23) 



ai = Ci + 



9l 

N 



1/N. Expanding (fTTj) - ()19|) to ffist order in l/rrib 

C2 CpT^Ois 



1a — C4 -|- 



1 + ^V. 

CpOts P-K,2 



+ 



A^ N 



'-H., 



7V7T.2 



rial = Tie, 



(24) 



(25) 



47r N 

We observe that to this order in the double expansion, the uncalculable power 
correction ~ -^^7777,3 does not appear in oi, to which it only contributes at order 
l/mfyN. On the other hand, the leading-twist effect ~ HTTn,2 is retained as well 
as the vertex corrections ~ V^^r- Both contribute to ai at order 1/A^. For a^ 
the hard-spectator term is absent altogether because it scales as 1/N^, but the 
nontrivial penguin loop corrections still contribute at order 1/A^. Since r'J^ag is 
already ~ l/rrib we omit all 1/A^ effects. 

We next show that within our approximation the expressions in ()24|) . (f^3j) 
receive no further corrections from weak annihilation. From jl] we recall that the 
annihilation coefficients bi appearing in (jllll can be written as 



61 



Cf 

ATS 



CiA\ 



Cf 



C^A\ + CeAi 



CsA\ + C5iAl + 4) + NCeA{ 



(26) 



(27) 



The parameters A^fl are not calculable in QCD factorization. However, they have 
been estimated in ^ from a diagrammatic analysis of the annihilation topologies 
in a way that keeps track of the correct counting in l/rub and 1/A^. One finds |3] 



AU ~ 1/A^ 



Aii^ 



l/nibN 



ta ~ N/mb 



f28) 



and therefore 

r^6i,4 ~ l/rUkN ta^ ~ l/ml (29) 

All annihilation effects thus contribute to r in (jlip only at second order in the 
double expansion, as anticipated above. Numerically the largest impact on r 
comes from 63, which in turn is dominated by the term ~ Cq. The contribution 
to 63 from C3 is highly suppressed ~ l/mi,N'^, the one from C5 even by l/mlN"^, 
and both are very small numerically. We remark that, unlike all other corrections 
to naive factorization, the annihilation term VAbs is leading in the large- A^ limit. 

We conclude that plj) and (J25|l give indeed the amplitude coefficients complete 
through first order in l/mi, and 1/A^. We stress again that in this approximation 
these quantities are fully calculable. In other words, the problematic corrections, 
uncalculable in QCD factorization, from higher-twist spectator interactions and 
weak annihilation are at least doubly suppressed in the combined heavy-quark, 
large- iV expansion. This observation shows that the large- A^ limit yields a useful 
organizing principle complementary to the l/rrif, expansion. If the large- A^ limit 
is not entirely unrealistic, the double expansion will provide an additional tool 
to improve our theoretical control over two-body hadronic S-decay amplitudes. 
Experience with similar decays of K and D mesons suggests that considerations 
based on large-A^ arguments can be a reasonable approach to these problems 
|2ti|27|28j . After all, as discussed in this paper, for some applications approximate 
results are sufficient and precise calculations are not necessarily required. 

It is interesting to evaluate the approximations (j23), ()25|) numerically and to 
compare with the full NLO QCD factorization results that include the estimates of 
uncalculable power corrections. Using default input parameters we quote three 
values for the various coefficients: The first, second and third numbers give, 
respectively, the result for QCD factorization, for the l/ruh, 1/N approximation 
in (j23)) ^25|1 ■ and for naive factorization. -l- 



ai = 


1.00 + 0.02Z 


1.02 + 0.02Z 


1.03 


al = 


-0.033 - 0.007Z 


-0.038 - 0.006Z 


-0.027 


■>l = 


-0.056 - 0.007Z 


-0.041 


-0.038 


r = 


0.107 


0.084 


0.068 


6 = 


0.150 


0.065 






(30) 



For comparison, the default value for the annihilation correction to the penguin 
amplitude is VAip^ + 264) = —0.010 in the model of |1]. We see that to first 
order in the double expansion we recover the QCD factorization result to a large 
extent. Since we use the full NLO coefficient C4 in a^, there is strictly speaking 
an ambiguity whether or not to include in al the C3 terms, which cancel a small 
part of the NLO scheme dependence. If we include C3 in the vertex and penguin 
correction part of a%, still neglecting hard spectator scattering, we find —0.036 — 

"'"The coefficients in naive factorization are ai = Ci + Ci_/ i)i/3, with leading-log values for 
the Ci, and bi = 0. 



0.007i instead of the -0.038 - 0.006i above, and (r, 0) becomes (0.081,0.083) 
instead of (0.084, 0.065). Apparently the 1/N approximation proposed here does 
not seem to be entirely unrealistic, at least within the standard QCD factorization 
framework at next-to-leading order. 

4 Unitarity Triangle from S and sin 2/3 

4.1 Determining p and f] 

In this section we discuss the determination of the unitariy triangle by combin- 
ing the information from S with the value of sin 2/3, which is known with high 
precision from CP violation measurements in S — ;> J /'^Ks- As we shall see, this 
method allows for a particularly transparent analysis of the various uncertainties. 
Both p and fj can be obtained, which fixes the unitarity triangle. A comparison 
with other determinations then provides us with a test of the standard model. 
The angle (5 of the unitarity triangle is given by 

cot/? = ""'^ (31) 

1 - sin^ 2j3 



The current world average 

sin 2/3 = 0.739 ±0.048 (32) 

implies 

r = 2.26 ±0.22 (33) 

Given a value of r, p is related to f/ by 

p = l-Tfi (34) 

The parameter p may thus be eliminated from S" in (0), which can be solved for 
f] to yield 



(1 ± tS){1 ± r cos 0) - ^(1 - 52)(1 + r2 + 2r cos 0) - (1 ± r^) V^ sin' (j) 

"^^ (l±r2)5 

(35) 

So far, no approximations have been made and eqs. (jH^ . (jH^ are still com- 
pletely general. The two observables r (or sin 2/?) and S determine f] and p once 
the theoretical penguin parameters r and (f) are provided. It is at this point 
that some theoretical input is necessary. We will now consider the impact of the 
parameters r and 0, and of their uncertainties, on the analysis. 

We first would like to point out that the sensitivity of fj in (jHSI) on the strong 
phase (p is rather mild. In fact, the dependence on enters in (jHK|l only at second 
order. Expanding in we obtain to lowest order 



l + rS-^l-S^ , , , , 
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Figure 1: CKM phase fj as a function of the mixing-induced CP asymmetry S in 
B — > vr+TT" within the standard model for sin 2/3 = 0.739. The dark (light) hand 
reflects the theoretical uncertainty in the penguin phase = 0.15 ± 0.25 (penguin 
amplitude r = 0.107 ± 0.031^. 



This result is corrected at second order in through 



Af] 



'l-^^ + r(l + r^)- 



[1 + tS) 



rcl? 



2(l + r2)5 



(37) 



This feature is very welcome since it is in particular the strong phase that is dif- 
ficult to calculate with good precision. Nevertheless, we know from factorization 
in the heavy-quark limit that the strong phase is suppressed either by a^, if it 
arises from hard scattering, or by t^QCDJ'^^h for soft corrections. This means that 
even if is not accurately known, it will have very little impact on f] as long as 
it is of moderate size. Since r is also small, the second order effect from in ()37j) 
is even further reduced. As an example, for S* = one finds 



Ar/ 



1 



r 



r 1 + r2 2 



(3^ 



For typical values r ^ 0.1, this implies that |Af/| < 0.01 for up to 45°, which is 
already a large phase. For < 20°, which is more realistic, one has a negligible 
shift |Af/| < 0.002. Consequently, the relation in ()36|) is most likely a very good 
approximation to the exact result. Note that apart from neglecting the phase 0, 
no approximation is made in p6|l . The resulting expression is strikingly simple. 
f] is essentially determined by the CP violating observables S and r. The only 
dependence on the penguin parameter r is through an overall factor of (1 + r). 
Typically r ?^ 0.1, as predicted in QCD factorization but indicated also by other 
approaches. The effect is again a fairly small correction. A 100% uncertainty on 
this estimate of r would translate into a 10% uncertainty in f\. 



10 



The determination of f/ as a function of S is shown in Fig. ^ which dis- 
plays the theoretical uncertainty from the penguin parameters r and (p in QCD 
factorization. 

In the determination of fj and p described here discrete ambiguities do in 
principle arise. One source is the well-known ambiguity in relating sin 2/? to a 
value of P, or equivalently r = cot/3. Apart from the solution shown in (jnU, a 
second solution exists with the sign of the square root reversed. It corresponds to 
a larger value of (3, incompatible with the standard fit of the unitarity triangle. 
An additional ambiguity comes from the second solution for f/, which is the result 
given in ()H5j) with a positive sign in front of the square root. This case may be 
considered separately, but will usually also yield solutions in conflict with other 
information on the CKM phases. 

4.2 Standard Model Prediction for S 

In this section we shall use theoretical information on r and based on QCD 
factorization to compute the value of 5* expected within the standard model. 
Using ()34|). one can write (0) in the form 



r/-r 1 + 



Q r, 1+r COS cf> I \ {l+rcos4>)~ / ("iQ^ 

^^ \l+rcos<f>'l 'J ^ IJ- -^ / j(i+^cog^)2 

Since the terms ~ r^ sin^ (p (< 10^'^) are very small, 5* is well approximated by 

1 + 2;^ 1 + r cos (p 



Taking 
and 



r = 2.26 ±0.22 r^ = 0.35 ± 0.04 (41) 

r = 0.107 ±0.031 (/) = 0.15 ±0.25 (42) 



we find from (j39j) 

5 =-0.59 tr,Ur) t'oitiv) ;H^(r) t',foi<t>) (43) 

We note that the uncertainty from the hadronic phase is negligible and the 
uncertainty from the penguin parameter r is rather moderate. The error from 
T or sin 2/5, which reflects the current experimental accuracy in this quantity, 
is considerably larger. The dominant uncertainty, however, is due to f], which 
for the purpose of predicting S has here been taken from a standard CKM fit. 
Clearly, the large sensitivity of 5 to f/ is equivalent to the fact that in turn fj has 
only a fairly weak dependence on S. This feature was already discussed in [50| . 
Eq. pnj) shows that the standard model prefers negative values for S, but it is 
difficult to obtain an accurate prediction. 

11 



4.3 Rigorous CKM Bounds with Minimal Hadronic Input 

We will now relax the constraints on the penguin parameters r and coming from 
direct theoretical calculations and derive rigorous bounds on p and fj without 
relying on any detailed information about hadronic quantities. Specifically, we 
shall only assume that the strong phase fulfills 

-^<*<2 (44) 

In view of the fact that is systematically suppressed in the heavy-quark limit, 
and that typically ~ 0.2 from QCD factorization, this assumption is very weak. 
As has been shown in [ITj, the following inequality can be derived from ()35p 
for - sin 2/? < ^ < 1 



V> (i + T^)S (l + '^cos^) (45) 

This bound is still exact and requires no information on the phase 0. (The only 
condition is that (1 + r cos0) is positive, which is no restriction in practice.) 
Assuming now (jl^ , we have 1 + r cos > 1 and 



r/>i±^f-^^I^ if -sm2f3<S<l (46) 

(1 + T^)S 

We emphasize that this lower bound on f/ depends only on the observables r 
and S and is essentially free of hadronic uncertainties. It holds in the standard 
model and it is effective under the condition that S will eventually be measured 
in the interval [—sin 2/3, 1]. Since both r and are expected to be quite small, 
we anticipate that the lower limit (jIUI) is a fairly strong bound, close to the 
actual value of f] itself (compare (jSHl))- We also note that the lower bound (jlUj) 
represents the solution for the unitarity triangle in the limit of vanishing penguin 
amplitude, r = 0. In other words, the model-independent bounds for fj and p 
are simply obtained by ignoring penguins and taking S = sin 2a when fixing the 
unitarity triangle from S and sin 2(3. 

Let us briefly comment on the second solution for f/, which has the minus 
sign in front of the square root in (jH^ replaced by a plus sign. For positive S 
this solution is always larger than (jSHI) and the bound (jlHI) is unaffected. For 
— sm2j3<S<0 the second solution gives a negative fj, which is excluded by 
independent information on the unitarity triangle (for instance from indirect CP 
violation in neutral kaons [ek))- The additional solution for f/ is illustrated in 
Fig. 121 for r = 0, the case relevant for the lower bound. 

Because we have fixed the angle (3, or r, the lower bound on fj is equivalent 
to an upper bound on p = 1 — rfj. The constraint ^^ may also be expressed as 
a lower bound on the angle 7 



7 > arctan ^ ' (47) 
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Figure 2: Discrete ambiguities for f] as a function of S with sin 2/5 = 0.739, r = 0. 
For — sin 2(3 < S < 1 the middle branch defines the lower bound on fj, which is 
not affected by the additional solution. 




Figure 3: Lower bound on rj as a function of S for various values of sin 2/? 
(increasing from bottom to top). 
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Figure 4- Lower bound on '-f as a function of S for various values of sin 2/5 
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Figure 5: Region in the (p, ff) plane constrained in a model-independent way by 
sin 2/5 = 0.739 ± 0.048 (shaded sector) and various possible values for S. The 
allowed area is the part of the shaded sector to the left of a given line defined 
by S . These lines correspond, from bottom to top, to 5=— 0.6, —0.3, 0, 0.3, 0.6 
and 0.9. The bound becomes stronger with increasing S. The result of a standard 
unitarity triangle fit (dotted ellipse, from [291) is overlaid for comparison. 
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or a lower bound on Rf 



1+tS- VT^^s^ 



Rt^^il- P? + V^ > ^^^;,^ (48) 

In Figs. El and ^ we represent the lower bound on f/ and 7 as a function of S 
for various values of sin 2(3. From Fig. IHl we observe that the lower bound on f/ 
becomes stronger as either S or sin 2(3 increase. The sensitivity to sin 2(3 is less 
pronounced for the bound on 7. Similarly to f], the minimum allowed value for 7 
increases with S. A lower limit 7 = 90° is reached for S = sin 2(3. 

In Fig. El we illustrate the region in the (p, fj) plane that can be constrained 
by the measurement of sin 2/5 and 5* using the bound in (jlF)|l . 

We finally note that the condition rcos0 > 0, which is crucial for the bound, 
could be independently checked [31j by measuring the mixing- induced CP-asym- 
metry in Bg — > K^K~ . This is because the hadronic physics of B^ — > K^K~ 
is related to B^ —>■ tt+tt^ by U-spin symmetry, a feature that has already been 
employed for CKM phenomenology [221 • Our purpose here is to use information 
from Bg -^ K^K~ in order to obtain additional input for B^ -^ tt+tt^ within 
the approach suggested above. To this end we write the CP-violation observable 
5", defined in analogy to (0), for the case of Bg — > K^K~ . Neglecting the small 
phase of Bs -Bs mixing, we find 



2?7 {kr cos (f) — p) 
{kr cos — PY + ff + k'^r^ sin^ 



siB.^K^K-)^ .:'"rr°i :".„,. m 



where 

1 — A^ 
k = ^^ ^ 20 (50) 

for Wolfenstein parameter A = 0.22. Because of the different CKM hierarchy of 
the 6 — > s transition, the penguin contribution ~ r is strongly enhanced by a 
factor of /c ~ 20 compared to the case of Bd —>■ tt^tt^. On the other hand, the 
purely hadronic quantities r and (p are identical to the corresponding parameters 
for Bd — > TT+TT" in the limit of exact U-spin symmetry. We notice that the 
presumably largest effects from U-spin breaking, coming from the difference of 
decay constants fx, fn and form factors for Bg —>■ K and Bd —>■ tt transitions, 
largely cancel in the ratio r of penguin over tree amplitudes. Explorations of 
further sources of U-spin breaking can be found in |31l33j . We shall assume that 
r «i 0.1 as indicated by QCD factorization. Then the relevant penguin parameter 
in ()49|) is kr ^ 2, which dominates the much smaller values for p ~ 0.15. As a 
consequence, (UHl) predicts, in the standard model, the sign of S{Bs — > K^K~) in 
correspondence with the sign of r cos 0. In QCD factorization this sign is positive 
and one expects 

S{Bs -^ K+K-) ^ r] (51) 

A future measurement of this observable will then provide a test of the assump- 
tion made to obtain the above bounds. We remark that also from other charmless 
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hadronic B decays a sizable penguin amplitude is required, independently of de- 
tailed QCD calculations, where r = 0.1 is a typical value. It is thus basically 
excluded that kr will be much below about 2 and that the term —p in the nu- 
merator of (j49j) will be able to compete so as to change the conclusion. On the 
other hand, an extreme value of the phase ~ 7r/2, for instance, could be in- 
dicated through the observation of a very small S{Bs -^ K~^K~), which would 
typically amount to a few percent. (In the approximation above one would ob- 
tain S{Bs -^ K^K^) K, —fip/2, but the mixing phase could then no longer be 
neglected. It would tend to further reduce the asymmetry.) 

In this comparison it is legitimate to assume the validity of the standard model 
throughout, as the strategy is to look for new physics via inconsistencies under 
this assumption, exploiting a multitude of experimental results and eliminating 
hadronic uncertainties. 



5 Unitarity Triangle from 5, Rh and sin2(/)(i 

In the previous section we have considered a determination of the unitarity trian- 
gle within the standard model. In particular we have assumed that the phase of 
Bd-Bd mixing /5, measured in the time-dependent CP asymmetry of i? — > J/'^Ks, 
is indeed the angle (3 of the CKM unitarity triangle, satisfying the relation 
tan/5 = fi/{l — p). In the presence of new physics this needs no longer to be 
the case. From this perspective a different analysis is of interest. To be spe- 
cific we shall assume the plausible scenario where the new physics contributions 
modify the phase of Bd-Bd mixing 0^, whereas the B decay amplitudes remain 
unchanged. The CP asymmetry in B —^ J /'^Ks ()H2j) must then be interpreted as 
the quantity sin20d. Since we can no longer relate sin20d to p and f/, we should 
fix it to the experimental value in (jH^ when using Q, where (3 is to be replaced 
by (pd- A similar analysis has already been carried out in H^ (see also j3]). This 
paper gives an interesting discussion of the possible implications of new physics 
for the case at hand, including a consideration of rare decay channels. The re- 
quired penguin-to-tree ratio is estimated in jTHj in a phenomenological way. Here 
we would like to rephrase the analysis in terms of the hadronic quantities r and 0, 
which are more directly related to theoretical input. We thus obtain formulas in 
line with our general approach to isolate CKM parameters and hadronic physics 
in a transparent way. 
Writing 



we have 



5" = Im 



p = Rb cos 7 ?7 = -Rfe sin 7 Rb = y p"^ + if (52) 

2is. {Rb COS •y + r cos (f) — iRb sin •y)^ + r^ sin^ ( 



{Rb cos y + r cos (py + Rl sin^ 7 -|- r^ sin^ 



(53) 



From this relation, for given values of r and 0, and using the experimental results 
for S, sin20d and Rb, 7 can be determined. Under the assumptions specified 

16 




Figure 6: Determination of the unitarity triangle from S, Rb and sm20rf. We 
take sin20rf = 0.739, Rb = 0.39, the penguin parameters r = 0.107 ± 0.031, 
(f) = 0.15 ± 0.25, and assume a hypothetical value of S = —0.4. The resulting 
ranges of the CKM angle 7 are 7 = 67° ± 4° (right band) and 7 = 138° ± 4° (left 
band) in this example and reflect the uncertainty in r and (p. The circular area 
shows the Rb constraint |^jj. 

above this will be the true angle 7 of the CKM matrix, independent of the B^-Bd 
mixing phase 0^. Once 7 is known, p and fj follow from (j52j) . 
Experimentally one has ;24^29j 



sm. 



0.739 ± 0.048 



Rb = 0.39 ± 0.04 



(54) 



As emphasized in [IB], there is a discrete ambiguity in the sign of cos2(j)fi, which 
yields two different solutions for 7. The larger value of 7 will be obtained for 
negative cos 20^. The analysis is represented in Fig. IHl assuming a particular 
scenario for illustration and displaying the impact of the theoretical uncertainty 
in r and (h. 



6 Direct CP Violation 

So far we have considered the implications of mixing- induced CP violation, de- 
scribed by S. In the following we shall investigate how useful additional informa- 
tion can be extracted from a measurement of the direct CP violation parameter 
C. An alternative discussion of this question can be found in jlOj. 

The observable C (see (0)) is an odd function of 0. It is therefore sufficient 
to restrict the discussion to positive values of 0. A positive phase is obtained 
by the perturbative estimate in QCD factorization, neglecting soft phases with 
power suppression. For positive also C will be positive, assuming r^ > 0, and a 
sign change in will simply flip the sign of C. 

In contrast to the case of S, the hadronic quantities r and play a prominent 
role for C, as can be seen in (jH)). This will in general complicate the interpretation 
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of an experimental result for C. One aspect of this can be seen as follows. It 
is usually expected that small values of the weak phase fj and the strong phase 
correspond to a small CP asymmetry C. However, in principle this need not 
be the case. As a counterexample let us consider the scenario where p = —0.35, 
fj = 0.07, r = 0.35 and = 0.2 ?» 11°. Of course this implies a very large angle 7, 
but this would be possible if the presence of new physics invalidates the standard 
unitarity triangle analysis (still the constraint from Ri, is obeyed). Although both 
f] and the strong phase are very small, these numbers give C = 0.995. More 
generally, such a situation occurs if p = — r, leading to a cancellation in the 
denominator of C. In this case, assuming also that is small, we get 

7]^ + [rcpY 

which takes its maximal value C = 1 for fj = rcf). Clearly, a scenario of this type 
requires a very peculiar coincidence and may seem unlikely. Nevertheless the 
example illustrates that the proper interpretation of C can be rather involved. 

The analysis of C becomes more transparent if we fix the weak parameters and 
study the impact of r and 0. An important application is a test of the standard 
model, obtained by taking p and f/ from a standard model fit and comparing the 
experimental result for C with the theoretical expression as a function of r and 
0. 

Let us first derive a few general results. An important question is the max- 
imum value of C, for given p and r], allowing an arbitrary variation of r and 
0. 

Varying r we find that C takes its maximum for 



r = R^ = ^p^ + fj^ (56) 

independently of 0. The resulting maximum Cmax{(p) at r = i?^ can be written 
as 

C^M = J'""""^^ (57) 

1 + COS 7 cos 

and only depends on and 7. Viewed as a function of it can reach its absolute 
maximum C = 1 for cos0 = — cos 7. 

A useful representation is obtained by plotting contours of constant C in the 
(r, 0)-plane, for given values of p and fj. This is illustrated in Fig. [7| for the 
standard model best-fit result p = 0.20, fj = 0.35 [29j. Within the standard 
model this illustrates the correlations between the hadronic penguin parameters 
r and and direct CP violation in i? ^ tt+tt^ decays. Upper limits on r and 
imply an upper limit on C unless they acquire unreasonably large values. For 
example, r < 0.15 and < 0.5 yield C < 0.21. 

We may relax the assumption of the validity of the standard model and dis- 
cuss the parameter C from a different perspective. Similarly to the previous 
section, we consider the rather general scenario where new physics renders the 
standard unitarity triangle fit to determine 7 invalid, while the extraction of Rb 
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Figure 7: Contours of constant C in the (r, ipj-plane for fixed p = 0.20 and 
f] = 0.35. 

and the B —>■ tt+tt" amplitudes remain essentially unaffected. In this situation it 
is convenient to slightly rewrite (jHI) as 



C 



2k sin 7 sini 



1 + K^ + 2k cos 7 cos( 



(58) 



where we have introduced n = r/Rb = \P/T\. If we treat 7 as unconstrained, we 
can still place an upper bound on C by maximizing C with respect to 7. Denoting 
this maximum by C we find 



C 



2k sin I 



1 + k'^Y ~ 4fi;2 cos^ I 



(59) 



where the maximum occurs at cos 7 = —2k cos 0/(1 + k^). 

If K = 1, or equivalently r = Rh, then = 1 independent of 0, and no useful 
upper bound is obtained. On the other hand, if k < 1, then C is maximized 
for = 7r/2. Under the general assumptions stated above and without any 
assumption on the strong phase we thus find the general bound 



C < 



2k 



1 + k2 



(60) 



For the conservative bound r < 0.15, k < 0.38 this implies C < 0.66. The 
bound on C can be strengthened by using information on 0, as well as on k, and 
employing (J3IJ|) . Then k < 0.38 and < 0.5 gives C < 0.39. 
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7 Tests of Factorization Predictions 

The analyses described above require theoretical input on the penguin parameter 
rexp(i0). We have relied on model-independent calculations based on the heavy- 
quark limit of QCD, including model estimates of subleading effects. To reinforce 
the validity of the approximations, in particular the large-m;, limit for realistic 
values of the 6-quark mass, it is important to test other predictions, obtained 
within the same framework, against experiment. For this purpose it is necessary 
to keep in mind that both hadronic effects as well as effects from new physics could 
in general be the origin of any discrepancy. Both effects need to be disentangled 
as far as possible. Especially useful tests of the QCD aspects in hadronic B decays 
are those that have no, or very little, dependence on weak phases and potential 
new physics contributions. We shall discuss several such tests, which pertain to 
the essential ingredients for the P/T ratio rexp{i(f)), namely the tree amplitude, 
the penguin amplitude and annihilation effects in S — > vrvr decays. 

7.1 Tree amplitude 

The first example is a factorization test for the rate of the tree-type decay B^ — >■ 
7r^7r° by taking a suitable ratio with the semileptonic rate of -B ^ nlu. For similar 
decays such observables have already been discussed e.g. in i34i. More recently 
it has been proposed to employ the decay B -^ -kIv for an estimate of the tree 
amplitude in B ^ tt^tt", assuming factorization [35j. This analysis is similar to 
the one suggested here, however our main emphasis is somewhat different. We 
consider the factorization test as a cross-check of dynamical calculations based on 
the heavy-quark limit (4j , rather than a way to determine the tree amplitude. The 
reason for this distinction is the fact that the tree contribution in B^ -^ -k^-k^ 
is not exactly the same as the one required for B -^ tt+tt^. A related discussion 
can also be found in [S^ . 

The differential decay rate for B^ -^ -k^I^v is given by 



ds 






192773 



ub\ 



>^ns)fw 



(61) 



with 



K{s) 



l + rl + s"^ -2s-2r^~ 2r^s 



niz 



m 



B 



m 



B 



(62) 



Here g^ is the invariant mass of the lepton pair and f+{q'^) is a. B ^ n transition 
form factor. Eq. ()61|) is valid for leptons I = e, fi where the mass is negligible. The 
pion mass effect is also very small, r^r = 7 x 10"^, and can likewise be neglected. 
In this case the branching fractions of B^ —>■ tt^tt" and B^ -^ ir'l'^v are related 
through 



B{B^ 



tt+ttO 



)=37r7^|K 



ud\ 



dB{B^ 



TT-I+U) 



dq^ 



t{B- 



r{B, 



\0'l+ 0.2] 



(63) 
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where ai, a2 are QCD coefficients |2p4| . In (jU^ the dependence on \Vub\ and /+(0), 
which are not known precisely, have cancelled out. Once the branching fractions 
are measured, |ai + 02] can be extracted experimentally via ()63|) and compared 
with the theoretical prediction. This test is useful since neither QCD penguins 
nor weak annihilation corrections affect the B~^ — * 7r"'"7r° amplitude. It can thus 
provide us with a check on the tree amplitude ai + 02, which includes non-trivial 
hard-spectator interactions in QCD factorization, but does not depend on the 
other complications. Moreover, this test is independent of CKM phases and very 
unlikely to be modified by non-standard physics. It probes, as desired, crucial 
aspects of the QCD dynamics in 5 — ;> vrvr decays. 

In order to determine the differential semileptonic branching ratio dBsi/dQ^ 
= dB{Bd — > 7i~ l~^ u) / dq"^ at g^ = 0, one needs to fit the q^ spectrum of the 
semileptonic decay. We may use the expression 



dBsL dBsL 



dq^ dq^ 



! + «:„ . 

2 (64) 



g2 



which follows from (jfjlll and the parameterization of the form factor suggested in 



/+(0) 1 



(65) 






At present the data on B^ -^ -n^l^v decays are not yet accurate enough to give a 
stringent test [33]. The situation should improve substantially in the future and 
will then yield valuable information on QCD dynamics m B ^ -kti decays. 

7.2 Penguin-to-Tree Ratio 

The decay mode B^ -^ tt^K^ is essentially a pure penguin process, up to a 
negligible rescattering contribution |3]. The ratio of the B^ -^ n^K^ to the 
B~^ — »• 7r"'"7r° branching fraction is therefore a useful probe of the pegnuin-to-tree 
ratio |4l36j . In analogy to the relevant parameter r in B —>■ tt+tt", one may define 
a quantity f, which can be expressed through observables: 



al + r^al + rAh)-KK 



ai + 02 



K 



ub 



Vr_ 



cb 



f. 



B{B+ -^ TT + ifO) 



fK\2B{B+^n+nO) 



0.099±0.011±0.008 



(66) 
Here CP-averaged branching fractions are understood. The quoted number is de- 
rived from current experimenatl results, where the ffist error comes from \Vub/Vcb\ 
|29j . the second from the branching ratios (see Table Ej). The parameter f differs 
from r in the numerator through a different annihilation correction (63 instead of 
63 -|- 264 (fTT|) ) and through small SU{3) breaking differences in the light-meson 
distribution amplitudes {ttK instead of tttt). In the denominator f has the pure 
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5° ^ TT + TT- 


B+ -^ tt+ttO 


5° ^ vrOvrO 


Acp(7r+7rO) 


4.6 ±0.4 


5.3 ±0.8 


1.90 ±0.47 


-0.07 ±0.14 


BO -^ K+TT- 


B+ ^ i^+7rO 


5+ ^ K%+ 


5° ^ i^o^o 


18.1 ±0.8 


12.8 ±1.1 


19.6 ±1.5 


11.2 ±1.4 


AcpiK+n-) 


Acp(i^+7rO) 


AcpiK^^) 


AcpiK'n^) 


-0.09 ±0.03 


0.00 ±0.07 


-0.01 ±0.06 


±0.03 ±0.37 


5° -^ K+K- 


fi+ ^ K+i^o 


BO ^ i^o^o 




<0.6 


< 1.3 


< 1.6 





Table 3: Current world average values for B —>■ Tm, Kn branching ratios (CP 
averaged, in units of 10^^) and direct CP asymmetries Acp{f) = (r(i? ^ /) — 
TiB^f))/iTiB^f) + TiB^f)) 



tree amplitude oi ± a2, while r has oi corrected by small penguin and annihila- 
tion terms. Despite these differences in details, the structure of r and r are very 
similar. In fact, the theoretical value for r from QCD factorization 



0.081 ± 0.016 ± 0.016 = 0.081 ± 0.023 



(67) 



is very close to the corresponding value for r, and both are in agreement with the 
experimental number in (jUBj) . 

A final comment concerns the branching ratio for B -^ ti^tt^, which appears to 
be somewhat larger experimentally (Table El) than expected in recent theoretical 
calculations |22j, even though the error bar is still large. It should be stressed that 
B — >■ 7r*^7r°, being colour-suppressed (amplitude involving 02), is highly sensitive 
to the dynamics of hard spectator interactions, which so far are only known to 
lowest order in QCD and depend on poorly known input within the factorization 
framework. These uncertainties strongly affect B -^ vr^vr^, but are considerably 
smaller in B ^ vr+vr'^ and B -^ vr+vr", as already pointed out in .2\. In [36j a 
scenario with an enhanced B -^ tt^tt^ rate, without the need for very unusual 
hadronic input, was suggested. Such a scenario with large 02 could be checked 
using the factorization test discussed in the preceding subsection. We emphasize, 
however, that the uncertainties in 02 specific to i? ^ tt^tt^ have essentially no 
impact on the penguin-to-tree ratio r exp(i0), because the dominant hadronic 
physics is characteristically different. Even a relatively large value oi B[B — »• 
Tr^Tr*^) does therefore not invalidate the theoretical results for r exp(z0). 



7.3 Annihilation Decays 

Amplitudes from weak annihilation represent power suppressed corrections, which 
are uncalculable in QCD factorization and so far need to be estimated relying on 
models Q. At present there are no indications that annihilation terms would 
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be anomalously large, but they do contribute to the theoretical uncertainty. Ef- 
fectively, annihilation corrections may be considered as part of the penguin am- 
plitudes. To some extent, therefore, they are tested with help of the quantity f 
discussed in the previous subsection. Nevertheless, in order to disentangle their 
impact from other effects it is of great interest to test annihilation separately. 
This can be done with decay modes that proceed through annihilation or at least 
have a dominant annihilation component. 

An example is the pure annihilation channel B^ -^ D~K^ . Even though this 
case is somewhat different from the reactions of primary interest here, because of 
the charmed meson in the final state, it is still useful to cross-check the typical 
size of annihilation expected in model calculations. Treating the D meson in 
the model estimate for annihilation jlj as suggested in 3], one finds a central 
(CP-averaged) branching ratio of B{Ba -^ DjK^) = 1.2 x 10^^ and an upper 
limit of 5 X 10~^[3HI. This is in agreement with the current experimental result 
(3.8 ± 1.1) X 10-5 (see refs. in [HH]). 

Additional tests should come from annihilation decays into two light mesons, 
such as the B -^ KK modes in Table El These, however, are CKM suppressed 
and only upper limits are known at present. The K^K^ and K^K^ channels have 
both annihilation and penguin contributions. On the other hand B —>■ K^K~ 
is a pure weak annihilation process and therefore especially important. Further 
discussions can be found in |4p36j . 

8 Conclusions 

In this paper we have proposed strategies to extract information on weak phases 
from CP violation observables in B ^ ti^ti^ decays even in the presence of 
hadronic contributions related to penguin amplitudes. Our main results can be 
summarized as follows. 

• An efficient use of mixing-induced CP violation in i? — > tt+tt^ decays, mea- 
sured by S, can be made by combining it with the corresponding observable 
from B -^ V'-ft's) sin 2/5 or r = cot/3. 

— The unitarity triangle parameters p and f/ can then be obtained in 
closed form as functions of the observables r, 5* and the hadronic 
penguin parameters r, (see eqs. (jSll), (ISSI))- 

— The sensitivity on the hadronic quantities, which have typical values 
r ~ 0.1, ~ 0.2, is very weak. In particular, there are no first-order 
corrections in 0. For moderate values of its effect is negligible. 

— Neglecting 0, the dependence of f/ on r comes merely through an overall 
factor {1+r). The impact of the uncertainty in r f» 0.1 becomes clearly 
visible and is seen to be greatly reduced. A simple determination of 
the unitarity triangle from r and 5* is thus possible (see eq. (j36|) ). 
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• The parameters r^, 1 — p, Rt and 7 are bounded from below, depending only 
on r and S and essentially without relying on hadronic input (see eq. (UHl))- 

• The parameter of direct CP violation C depends much stronger on hadronic 
input than S", but yields complementary information and can constrain r 
and within the standard model. 

• An alternative analysis of S and C especially suitable for the presence of a 
new-physics phase in B^-Bd mixing was discussed. 

As an input to the phenomenological discussion we also studied the calculation 
of the penguin parameters r and cj) in QCD: 

• We have analyzed r and within QCD factorization with a particular view 
on theoretical uncertainties. 

• B —>■ TT+TT^ amplitudes can be expanded simultaneously in l/rub and 1/A^, 
which leads to an interesting pattern of simplifications. All power cor- 
rections suffering from infrared endpoint divergences in the factorization 
formalism are at least of second order in this double expansion, while the 
most important effects survive at linear oder in l/rrif, or 1/A^. 

• The different contributions to -B ^ tt+tt^ amplitudes, the tree-component, 
the penguin-to-tree ratio and annihilation effects, appear in similar form 
in other B decay channels, such as B^ — »• tt+tt", B"^ —>■ tt^K^ and B — > 
K~^K~. These can be used to validate theoretical predictions, separately 
for the various components. 

The results presented in this paper should be useful for interpreting the forth- 
coming experimental measurements of CP violation in 5 — > 7r"'"7r~ decays in a 
transparent way and help to achieve a reliable control over theoretical uncertain- 
ties. 
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